Universality of anisotropic fluctuations from numerical simulations of turbulent flows 
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We present new results from a direct numerical simulation of a three dimensional homogeneous 
Rayleigh-Benard system (HRB), i.e. a convective cell with an imposed linear mean temperature 
profile along the vertical direction. We measure the SO(3)-decomposition of both velocity structure 
functions and buoyancy terms. We give a dimensional prediction for the values of the anisotropic 
scaling exponents in this Rayleigh-Benard systems. Measured scaling does not follow dimensional 
estimate, while a better agreement can be found with the anisotropic scaling of a different system, 
the random-Kolmogorov-flow (RKF) [1] . Our findings support the conclusion that scaling properties 
of anisotropic fluctuations are universal, i.e. independent of the forcing mechanism sustaining the 
turbulent flow. 



Small scales turbulent statistics is a challenging open 
problem for both theoretical and experimental studies 
in hydrodynamical systems [2]. Typical questions are 
connected to the understanding of the universality issue, 
i.e. to which extent small-scale turbulent fluctuations are 
statistically independent of the large-scale set-up used 
to inject energy in the flow. Robustness of small-scale 
physics cannot be exact. For instance, different forcings 
may inject large-scale different anisotropic fluctuations, 
which must have some direct/indirect influence on small- 
scale statistics. 

A first strong requirement for universality to hold 
is therefore that large-scale anisotropic fluctuations be- 
comes more and more sub-leading by going to smaller and 
smaller scales. In other words, at scales small enough, the 
omnipresent and universal isotropic fluctuations must be 
the leading statistical components. Such a requirement is 
always observed in both experiments and numerical simu- 
lations, although some subtle effects may show up due to 
the existence of anomalous anisotropic scaling (see [1,3- 
5] for a detailed discussion of this issue). Another im- 
portant question which must be asked about universality 
of small scales statistics, is connected to the anisotropic 
components on their own, independently on their com- 
parison with the isotropic ones. In particular, it is impor- 
tant to understand whether the anisotropic components 
of any turbulent correlation functions have a scaling be- 
havior characterized by universal exponents or not, in 
the limit of high Reynolds numbers. 

In this Letter we present first results of an attempt to 
study the small-scale anisotropic behavior of a homoge- 
neous three dimensional Rayleigh-Benard system (HRB), 
i.e. a convective cell with fixed linear mean temperature 
profile along the vertical direction. The main focus of our 
analysis is a comparison between the statistical behavior 
of HRB system with a completely different anisotropic 
flow, a random-Kolmogorov-flow (RKF) [1,3]. From the 
comparison, we show that the two systems have almost 
indistinguishable, in the limit of our numerical resolu- 



tion, small-scale anisotropic (and isotropic) scalings, i.e. 
we find a high degree of small-scale universality for all 
measurable anisotropic components. This result is par- 
ticularly relevant because its validity is only possible if 
HRB has anomalous (to be defined below in details) 
anisotropic small-scale fluctuations. 

This Letter is organized as follows. First we briefly 
discuss the physics of HRB flow and the details of our nu- 
merical simulations. Second, we review the technique of 
SO (3) decomposition to disentangle different anisotropic 
contributions to any velocity correlation functions. We 
then present our numerical results on the HRB. 

We first show that the observed anisotropic scaling is 
anomalous, i.e. it does not follow the dimensional predic- 
tions than can be derived by an analysis of the equation 
of motion. We then address in details small-scale univer- 
sality by making the comparison between HRB and RKF 
anisotropic properties, the central point of the present 
Letter. 




FIG. 1. A snapshot of the temperature fluctuations in the 
3D HRB system. The system is fully periodic. Notice the 
presence of typical convective structures [6] as the neat plume 
on the bottom/left of the picture. 
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An Homogeneous Rayleigh-Benard system (HRB) is a 
convective cell with fixed linear mean temperature pro- 
file along the vertical direction. The flow is obtained 
by decomposing the temperature field as the sum of 
a linear profile plus a fluctuating part, T(x,y, z;t) = 
T'(x,y,z;t) + (AT/2 - zAT/H), where H is the cell 
height and AT the background temperature difference. 
The evolution of a HRB system can be described by a 
modified version [14] of the Boussinesq system [6]: 

d t v + (v ■ V)v = -Vp + vV 2 v + agT'z (1) 

AT 

d t r + (v ■ v) r = K v 2 r - — Vz . (2) 

H 

Fully periodic boundary conditions are used for the ve- 
locity field, v, and temperature, T', fields. 
For large Rayleigh numbers, HRB shows a turbulent con- 
vective dynamics with absence of both viscous and ther- 
mal boundary layers [14]. The Bolgiano scale, Tb = 
e 5/4 N~ 3/4 (agy 3/2 , is of the order of the integral scale of 
the cell (H), hence temperature fluctuations have a lead- 
ing role only at the largest scales in the system. This 
has been already shown in a similar simulation [8], and 
is consistent with the picture presented in [7]. 

The main advantage of the HRB system is that the 
intrinsic homogeneity along the three directions allow a 
systematic study of scaling properties without spurious 
(non-homogeneous) effects, always present in standard 
Rayleigh-Benard systems with boundary layers. 

In order to asses the importance and properties of 
anisotropic components for any correlation function it is 
necessary to make a decomposition on the complete ba- 
sis of the SO (3) group [9]. In particular, in the following, 
we are mainly interested in the SO (3) decomposition of 
scalar quantities, as for the case of velocity longitudinal 
structure functions, S( p \r) — ( [(v (r) — v (0)) • r] p ) : 

oo j 

S {p) (r)=Y, E S ^(r)Y jm (f); (3) 

j=0 m=—j 

where the indices (j, m) label the total angular momen- 
tum and its projection on a reference axis of the spherical 
harmonics Yj m (f), respectively (see [3,9] for more de- 
tails). The physics is hidden in the projections, S^(r). 
We are interested to measure what are the scaling prop- 
erties (if any) of each projections on different anisotropic 
sectors: 

Sfi{r)~c om r^v\ (4) 

where we have assumed, on the basis of theoretical [9] 
and numerical [3] evidences, that the scaling exponents 
do not depend on the m index. To go back to the univer- 
sality issue discussed at the beginning; we expect that the 
coefficients cj m are strongly dependent on the anisotropic 



properties of the large-scale physics while the values of 
the scaling exponents, £ J (p), must enjoy a much higher 
degree of universality. In other words, whether a partic- 
ular sector is alive, Cj- m | > 0, or not, Cj m — 0, depends 
on the forcing; while, once that sector is switched-on, the 
way it propagates to small-scale is forcing- independent. 
This picture can be proved on a rigorous basis for some 
problems of scalar/ vector advection by Gaussian, white- 
in-time, velocity fields (Kraichnan models [10]). 
Concerning the SO (3) analysis let us notice that veloc- 
ity structure functions have even parity with respect to 
r, therefore projections on odd j values vanish. In the 
following, we consider also mixed velocity and tempera- 
ture structure functions which have, on the other hand, 
dominant odd parity. 

From equation (1) one may easily write down the sta- 
tionary equation for the second order velocity structure 
functions; the extension of Karman-Howarth equation in 
the presence of a buoyancy term [11]. The result is, ne- 
glecting for simplicity tensorial symbols: 

(Sv(r) 3 ) ~ er + agzr ■ (ST(r) Sv(r)) (5) 

.3=0,1,... .7=0 j=l ® .7=1,2,... 

where with e we denote the energy dissipation and with, 
(8v(r) 3 ) and (ST(r) 8v(r)), the general third-order veloc- 
ity correlation and temperature-velocity correlation, re- 
spectively. The two terms on the r.h.s. of equation (5) are 
called respectively the energy-dissipation term and buoy- 
ancy term. In (5) we report for each term the value of its 
total angular momentum, j. Let us notice that the energy 
dissipation term in (5) has a non- vanishing limit, for high 
Reynolds numbers, only in the isotropic sector, j = 0. On 
the other hand, the buoyancy coupling, agz, brings only 
angular momentum j = 1. Due to the usual rule of com- 
position of angular momenta we have that the buoyancy 
term, agz ■ (ST(r) Sv(r)), has a total angular momen- 
tum given by the rule: j tot = 1 <S> j = {j - 1, j, j + 1}. 
Using the angular momenta summation rule for j we can 
decompose the previous equation obtaining the following 
dimensional matching, in the isotropic sector: 

(Sv(r) 3 ) 1=Q ~er + agzr (Sv(r) 5T{r)) J=1 + ... (6) 

and in the anisotropic sectors, j > 0: 

(Sv(r) 3 )^ agzr (6v(r) ST(r)) (j _ 1)+ ... (7) 

where only dominant contributions are reported. 

In the isotropic sector the buoyancy term is sub- 
dominant with respect to the dissipation term at scales 
smaller than the Bolgiano length, r < Lb [12]. Therefore, 
in our simulation the isotropic velocity fluctuations are 
closer to the typical Kolmogorov scaling, 5v(r) ~ r 1 / 3 , 
rather than to the Bolgiano-Obhukhov scaling, Sv(r) ~ 

r 3/5_ 

Let us now focus on the anisotropic sectors. Equation (7) 
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is the simplest dimensional prediction one can derive for 
this system consistently with the anisotropic properties 
of the buoyancy term, sector by sector. It plays a key 
role in the following because we will show that the ob- 
served anisotropic scaling in our HRB system differ from 
the matching (7), i.e. we measure anomalous anisotropic 
scaling exponents. 

Our HRB simulation was performed using a Lattice 
Boltzmann scheme, with spatial resolution of 240 3 . We 
stored roughly 270 statistical independent configurations. 
The Prandtl number for the simulation is equal to unit, 
and the Rayleigh number Ra — (agATH 3 )/(vK) = 
1.38 -10 7 . Measured Bol giano scale is Lb ~ 370, roughly 
one and half the cell size, while ag used in the equa- 
tion of motion (1) is 2 • 10~ 3 . A typical snapshot of the 
temperature field is shown in Figure 1. Notice the well 
detectable structures typical of all other Rayleigh-Benard 
cell [6,13,15,16]. In particular, there is a beautiful hot 
plume on the central bottom/left part of the picture. 

We now present our numerical results. In or- 
der to check the small-scale properties of the HRB 
system we have carried out the SO(3) decomposi- 
tion of both longitudinal velocity structure functions 
(3) up to order p — 6 and of the simplest set of 
buoyancy-like terms made of q velocity longitudinal in- 
crements and of one temperature increment, B( q,1 \r) = 
([(«(r)-,(0)).f]'(T(r)-T(0))>: 



B<«.U(r) = £ £ B^\r)Y jm (f). 



(8) 



j—0 m- 



The dimensional matching made in (7) can be extended 
to any order of correlation function, giving, in terms of 
the velocity and buoyancy projections, the leading scaling 
contribution: 



S (p) (r) ~ rB 



(P-2,1) 
j-l,m 



(r). 



(9) 



Denoting with x J (g, 1) the anisotropic scaling properties 



of the buoyancy terms, (r) ~ r^ 9 ' 1 ) we get the 

dimensional estimate linking velocity and buoyancy scal- 
ing: 



e» = i+xr 1 (p- 2 ' 1 ) 



(10) 



where we have added a subscript d to remind the reader 
that it is the result of a dimensional matching. 

Let us first discuss the issue of anisotropic anomalous 
scaling by making a comparison between the numerical 
measurements and the dimensional estimate (10). In Fig. 
2 we show a comparison between the velocity and buoy- 
ancy j = 4 projections for p = 3, 5. 
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FIG. 2. Log-log plot of quantities entering in the dimen- 
sional matching (9) for some anisotropic sectors vs r. Two 
top curves refer to Sj P g(r) (*) and to the buoyancy term 

rBfjf^ir) (□) with p = 3 and j = 4. The two bottom 
curves refer to the case p — 5 with symbols, (+) for the ve- 
locity projection and (x) for the buoyancy term. Inset: same 
for j = 6. Notice that the buoyancy term decays faster in 
all cases. In both plots, buoyancy terms are shifted along 
y-axis for the sake of presentation and we limited the r-range 
extension to those values with a statistically significant signal. 

The inertial scaling measured for the projection of the 
velocity structure function, <S (r) , is more intermittent 



than the corresponding buoyancy term, rB 



. (p-2,i) 

j — l,m 



(r). In 



other words, the observed velocity scaling is different 
from the dimensional matching derived from the equa- 
tions of motion: it is anomalous. This result holds for 
all measurable orders also for j = 2 and for j = 6 sector 
(j = 6 is shown in the inset). 

Let us now discuss the other important issue of uni- 
versality. We argue that not only HRB has anomalous 
anisotropic scaling but also that the measured behavior 
is indistinguishable from what measured in the random- 
Kolmogorov-flow [3]. The point is far from being triv- 
ial and must not be underestimated. The HRB has an 
anisotropic forcing, given by the buoyancy term, which 
acts at all scales, ~ gzST(r), i.e. there is also a di- 
rect inject of anisotropics at small-scale, at variance with 
the RKF where the forcing was only at large scales. In 
Fig. (3) we show an overall comparison of <5>j^(r) mea- 
sured on the HRB and on the RKF [3]. Comparison is 
limited to j = 4 and j = 6, because RKF data from 
the simulation of [3] have a sign inversion in the j = 2 
sector which makes comparison inconclusive. As can be 
seen the agreement is quite satisfactory, except for the 
very small scales, smaller then the viscous scale, where 
as usual the SO (3) decomposition suffers from interpola- 
tion errors. The small discrepancies at large scales are 
also to be expected: the inertial properties of the two 
flows have to match quite different conditions at large 
scale. 
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FIG. 3. Log-log plot of anisotropic j = 4 projections 
Sj P d( r ) vs r f° r both HRB and RKF flows. Each couple of 
-almost- coinciding curves refer to the comparison of HRB 
and RKF projections at a given order p — 3, 4, 5, 6 (top to 
bottom). Inset: same for j = 6. 

The same comparison for j = 6, shown in the inset, 
qualitatively supports the same result. 

The fact that incrtial-scales fluctuations of both flows 
are almost indistinguishable is the first important con- 
firmation of the universality of anisotropic fluctuations 
in sectors with j = 4, 6. Similar conclusions can be 
drawn for the j = 2 sector in different experimental set- 
up [5, 17, 18] (the only sector measurable, indirectly, in 
experiments) . 

Let us conclude by summarizing the two main results 
of this Letter. First, anisotropic fluctuations in Rayleigh- 
Benard systems are anomalous. Second, notwithstanding 
the direct influence of the forcing mechanism at all small- 
scale, anisotropic fluctuations are universal, i.e. the 
small-scale dynamics is dominated by anomalous fluctu- 
ations, coming from the self-organization of the inertial 
evolution. Similar behavior is at the origin of anoma- 
lous scaling in Kraichnan models of passive/vector ad- 
vection as already discussed in the introduction. In the 
latter case, one connects rigorously the anomalous iner- 
tial scaling with the existence of zero-modes of the inertial 
operator (see for example [19, 20] for a detailed analysis 
of anisotropic scaling in passive advection of scalar and 
vector quantities, respectively). Here, for Navier-Stokcs 
equation, one may only stress the striking similarities, 
without being able to push it to some rigorous statement. 
Concerning universalities for the isotropic scaling of this 
Rayleigh-Benard system, we notice that due to the large 
value of the Bolgiano scale -of the order of the system 
size- we expect to observe small-scale isotropic fluctua- 
tions close to the usual Kolmogorov 1941 scaling (plus 
intermittency). This is indeed the case. The Bolgiano- 
Obhukhov isotropic regime with 8v(r) ~ r 3 / 5 cannot be 
accessed from this simulation. In the framework of the 
dimensional matching (6) the existence of a Bolgiano- 



Obhukhov scaling in the isotropic sector corresponds to 
a leading influence of the buoyancy term in the scaling 
range [21]. 
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